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We develop an unambiguous and practical method to calculate one-loop 

quantum corrections to the energies of classical time-independent field config- 

^ ■ urations in renormalizable field theories. We show that the standard pertur- 

bative renormalization procedure suffices here as well. We apply our method 

to a simplified model where a charged scalar couples to a neutral "Higgs" 

field, and compare our results to the derivative expansion. 
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I. INTRODUCTION 



Solitons arise in many model field theories [T|] where nontrivial time-independent solu- 
tions to the classical field equations exist. The fluctuations of quantum fields about classical 
configurations are difficult to study and can qualitatively affect, even destabilize, solitons. 
Topological arguments support the stability of some particularly interesting solitons. How- 
ever, much less is known about the fate of "non-topological solitons" that appear as minima 
in the classical action, but have no deeper claim to stability when quantum effects are taken 
into account. 

The quantum corrections to the energies of classical field configurations are typically 
highly divergent in 3 + 1 dimensions. In some cases, like the Skyrme model 0, the un- 
derlying theory is non-renormalizable, so the quantum contribution to the soliton's energy 
is unavoidably cut-off dependent and ambiguous, like any other radiative correction in a 
non-renormalizable theory. 

In this Letter we describe a systematic and efficient procedure for calculating the quan- 
tum fluctuations about time- independent field configurations in renormalizable field theories. 
We show that all divergences can be removed by the same renormalization procedure that 
renders the perturbative sector of the theory finite. The only "ambiguities" are the well 
known scheme and scale dependences of the renormalization prescription that are resolved 
completely in the perturbative sector. The result of our program is a renormalized quantum 
"effective energy" , whose non-trivial minima (if they exist) describe solitons in the quantum 
theory. 

First we show how to regulate and renormalize the divergences in the sum over quantum 
fluctuations. Then we develop calculational methods that are efficient and practical enough 
that quantum effects can be included in a search for stable field configurations. Our results 
take the form of an effective energy, £(<ft(x), m, {g}), depending on the "profile function" of 
the renormalized field, <f)(x), the renormalized mass, m, and various renormalized couplings, 
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{g}, defined in the usual perturbative sector of the model. One may then search over the 
parameter space characterizing <p{x) for minima of £ while holding m and {g} fixed. Here 
we treat the simple case of a charged scalar coupled to the field 0. Our methods can be 
generalized straightforwardly to models including fermions, gauge fields and self-coupled 
scalars. However, our approach is limited in that we work only to order h, and we only 
consider spherically symmetric profile functions. 

The possibility that the top quark in the standard model might be described as a non- 
topological soliton [[J provided the original motivation for our work. The top quark Yukawa 
coupling, g t , leads to its mass, g t v, where v is the Higgs vacuum expectation value. For large 
g t it would therefore appear favorable to suppress the Higgs condensate in the vicinity of 
the t-quark, and the top quark would be a sort of "bag". However, gradient and potential 
energy terms in the Higgs sector of the classical action oppose the creation of a such a bag 
in the Higgs condensate. In order to study the problem at the quantum level it is necessary 
to regulate and renormalize the divergences in the i-quark fluctuations about a deformation 
of the Higgs condensate. It is crucial to hold the renormalized parameters of the standard 
model fixed while varying the possible profiles of the Higgs field. 

Bagger and Naculich studied this problem by making a derivative expansion Q. They 
also worked in a large- N approximation, in which there are no quantum corrections above 
order %. We would expect that if the t-quark is a bag, then the Higgs field will vary 
significantly — A0 ~ v — over distance scales of order the Compton wavelength of the top 
quark — A ~ 1/gtV- However, the t-quark mass, g t v, also sets the scale for the derivative 
expansion. Thus all derivatives are of the same size, making the expansion unreliable. Our 
method is designed for such situations. 

We can trace elements of our approach back to Schwinger's work on QED in strong fields 
||. Schwinger studied the energy of the electron's quantum fluctuations — the "Casimir 
energy" - - in the presence of a prescribed, static configuration of electromagnetic fields. 
He isolated the divergences in low orders of perturbation theory. Our work can be viewed 
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as an extension of Schwinger's to situations where the field is determined self-consistently 
by minimizing the total energy of the system including the Casimir energy. In addition we 
complete the renormalization program and develop practical computational methods in three 
dimensions. Dashen, Hasslacher and Neveu renormalized the divergent contributions to the 
energy of the 4 kink and sine- Gordon soliton in 1+1 dimensions using a simple version of the 
method we propose here |J. Ambiguities in these models recently pointed out and studied by 
Rebhan and van Nieuwenhuizen @ can also be resolved with our methods. Studies of solitons 
in renormalizable models often note that the divergences in the quantum contribution to the 
soliton energy can be cancelled by the available counterterms [[]]. However, we are not aware 
of any work in 3+1 dimensions in which renormalization of the field configuration energy 
is done in a manner consistent with on-shell mass and coupling constant renormalization in 
the perturbative sector. 



II. FORMALISM 

We consider a renormalizable field theory with a real scalar field coupled to a charged 
scalar ip. We take the classical potential V((f>) oc (0 2 — v 2 ) 2 , and if) acquires a mass through 
spontaneous symmetry breaking. At the quantum level we put aside the <fi self-couplings 
and consider only the effects of the <f) — ip interactions. We further restrict ourselves to 0(h) 
effects in the quantum theory, which correspond to one-loop diagrams. 

Our model is defined by the classical action 

S[(f),lf)} = J d 4 X |^(^0) 2 - _ y 2 )2 + q^q^ _ (1) 

+ a{d,<P) 2 -b{<p 2 -v 2 )-c{<p 2 -v 2 ) 2 } , 

where we have separated out the three counterterms necessary for renormalization and 
written them in a convenient form. At one-loop order in if), these are the only counterterms 
required. 

We quantize around the classical vacuum = v and define h = — v , so that 
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f 1 TfP 1 

S[h, if)) = I d 4 x\ -{d^hf -{h 2 + 2vh) 2 + d^d^ip - MVV> - #(^ 2 + 2vh)ip*ip (2) 

I 2 



+a{d^h) 2 - 6(/i 2 + 2/iv) - cih 2 + 2hv) 2 \ 



where M = ^fgv is the if> mass and m 2 = Xv 2 /3 is the h mass. 

The one-loop quantum effective action for h is obtained by integrating out if) to leading 
order in h. We are interested in time-independent field configurations h = h(x), for which 
the effective action yields an effective energy S[h] that has three parts: 

s[h) =e d [h] + e ct [h] + s^[h) , (3) 

where £ c \[h] is the classical energy of h, 
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£ cl [h) = J d 3 x U|V/f + ^(h 2 + 2vh) 2 j , (4) 

£ ct [h] is the counterterm contribution, 

S ct [h] = J d 3 x {a\Vh\ 2 + b{h 2 + 2hv) + c{h 2 + 2hv) 2 ) , (5) 

and S^[h] is the one-loop quantum contribution from if). S ct [h] and £$[h] are divergent, but 
we will see explicitly that these divergences cancel for any configuration h(x). 

We fix the counterterms by applying renormalization conditions in the perturbative sector 
of the theory. Having done so, we have defined the theory for all h(x). We choose the on-shell 
renormalization conditions 



2, _ n . rfS 2 



Si = 0, S 2 (m ) = 0, and 



0, (6) 



dp 2 

where Si and ^(p 2 ) are the one- and two-point functions arising only from the loop and 
counterterms as seen in Fig. |IJ. 
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FIGURES 
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(2gv) 2 U{p 2 ) gn 
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FIG. 1. One- loop diagrams. 

We denote the one-loop diagrams with one insertion by Q and with two insertions by 
U(p 2 ), and find 

Si = 2vgQ + 2vb , 
S 2 (p 2 ) = (2vg) 2 U(p 2 ) + gQ + b + (2v) 2 c + ap 2 . (7) 

Defining 

nV) ^ , (8) 

the renormalization conditions eq.(|6]) then yield 

a = -(2vg) 2 W(m 2 ) , b = -gSl , c = g 2 (m 2 W(m 2 ) - Ii(m 2 )) , (9) 

which we then substitute into the counterterm energy, eq. (^|). 

Now we consider the calculation of £^ [h] . This energy is the sum over zero point energies, 
2%uj, of the modes of if) in the presence of h(x), 

£M=H E M (io) 

a 

where E a are the positive square roots of the eigenvalues of a single particle Hamiltonian, 
H, given by 
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H = -V 2 + M 2 + g(h 2 + 2vh). (11) 

The fact that if) is complex accounts for the absence of \ in eq. (fL0|). 

is highly divergent. However our model is renormalizable and therefore the counter- 
terms fixed in the presence of the trivial h must cancel all divergences in £^. Rather than 
attempt to regulate the divergent sum in eq. (]T0|) directly, we study the density of states that 
defines the sum. We can isolate the terms that lead to divergences in £^ and renormalize 
them using conventional methods. 

For fixed h(x) the spectrum of H given in eq. (|TT|) consists of a finite number (possibly 
zero) of normalizable bound states and a continuum beginning at M 2 , parameterized by k, 
with E(k) = +y/k 2 + M 2 . Furthermore, H depends on h only through the combination 

X = h 2 + 2hv, (12) 

so we can consider £^ to be a functional of \. We restrict ourselves to spherically symmetric 
h. Then 

^M = E^ + E( 2 ^ + 1 ) Jdkp e (k)E(k) (13) 

3 I 

where pe(k) is the density of states in k in the £ th partial wave and the Ej are the bound 
state energies. pe(k) is finite, but the sum over £ and the integral over k are divergent. 
Furthermore 

Pi (k)= P r c (k) + l^i, (M) 

where 5e(k) is the usual scattering phase shift for the £ th partial wave, and pf ee (k) is the free 
(g = 0) density of states. This relationship between the density of states and the derivative 
of the phase shift is shown for example in ||. 

At the outset, we subtract p lree {k) from the density of states since we wish to compare 
£tp[x] to ^[0]. Viewing £^[x] as the sum of one loop diagrams, we see that only the diagrams 
with one or two insertions of g\ are divergent. A diagram with n insertions corresponds 
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to the n th term in the Born expansion, so all possible divergences can be eliminated by 
subtracting the first and second Born approximations from the phase shifts that determine 
the density of states. Standard methods allow us to construct the Born approximation for 
the phase shifts |§, which is a power series in the "potential" g\- 
We define the combination 



5i(k) = 6t(k) — 6jp(k) — 8f'{k) 



(15) 



where 8^{k) and S ( /'(k) are the first and second Born approximations to Si(k). We then 
have 



f(2), 



(ftp _ 



(27 



(16) 



+g 2 1 7^-M-P 2 )\m\ 2 



(2tt) ; 



where 



X(P) = / (ftxx{x)e 



-ip-x 



(17) 



and likewise for h(p). Both h and x are rea l an d depend only on q = \p\ for spherically 
symmetric h. We have subtracted out the order g and g 2 contributions by using 5t (k) instead 
of 5g(k), and added them back in by using their explicit diagrammatic representation in terms 
of the divergent constant Q and the divergent function U(p 2 ). 
We can now combine and Set and obtain a finite result: 

r°° 1 dS e (k) 

J 

where 



re 

^ + ^ct = E^ + E( 2 ^ + 1 ) / 

o J 



7T (i/c 



-E(k)+r 2 [ 



(18) 



g 2 <ig 
2^ 



(n(-g 2 ) - n(m 2 ) + m 2 n'(m 2 )) x(q) 2 + Av 2 q 2 W(-q 2 )h(qf 



(19) 



II is log divergent, but both {U(—q 2 ) — IT (m 2 )} and IT' are finite, so T 2 [h] is finite as well. 

Each term in the Born approximation to the phase shift goes to zero at k = 0, so by 
Levinson's theorem 5e(0) = 5e(0) = irni where ni is the number of bound states with angular 
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momentum £. As k — > oo, 5e(k) falls off like r, <^(fc) falls off like |, and 5^ (&) falls off like 
p. Since the Born approximation becomes exact at large k, 5e(k) falls like p-. Thus we see 
that the first subtraction renders each integral over k convergent. The second subtraction 
makes the i-sum convergent. We are then free to integrate by parts in (|l8|), obtaining 

S[h] = S cl [h] + T 2 [h] - - £(2£ + 1) / dk St(k)-=r-r + Y,{E 5 - M) . (20) 

7T g JO H/[K) j 

In this expression we see that each bound state contributes its binding energy, Ej — M, so 
that the energy varies smoothly as we strengthen h and bind more states. 

As noted in the Introduction, the representation of the Casimir energy as a regulated 
sum/integral over phase shifts plus a limited number of Feynman graphs was derived by 
Schwinger for the case of a prescribed background field. Our aim is to develop it into a 
practical tool to study the stability of non-trivial field configurations h(r). 



III. CALCULATIONAL METHODS 



In this Section we describe the method that allows us to construct £ [h] as a functional of 
h and search for stationary points. We now consider the calculation of each of the terms in 
eq. ( |2I1| ) in turn. The classical contribution to the action is evaluated directly by substitution 
into eq. ([|). T 2 [h] of eq. flTD| ) is obtained from a Feynman diagram calculation, 
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(47T) 



q 2 dq 
Jo 



dx< 



log 



M 2 + q 2 x{l - x) 



m 2 x{\ — x) 



M 2 — m 2 x(l — x) M 2 — m 2 x(l — x) 

X\l — x) A 2 2Zf \ 2 1 



m 2 (2i) 



Av 2 q 2 h(qY 



M 2 - m 2 x(l — x) 

The partial wave phase shifts and Born approximations are calculated as follows. The 
radial wave equation is 



- u", + 



+ 9x(r) 



U£= k Ui, 



(22) 



where k 2 > 0, and x{ r ) ~ > as r — ► oo. We introduce two linearly independent solutions to 
eq. (p2[), ul\r) and uf\r), defined by 



9 



u { l\r) =e l ^ k ' r) rh { l\kr) (23) 
u f\r) = e - ip t^rhf\kr) 

where hg is the spherical Hankel function asymptotic to e lkr jr as r — > oo, hf\kr) = 
(kr), and r) — > as r — > oo, so that (r) — > e lfcr and u^(r) — > e~ lfcr as r — > oo. 
The scattering solution is then 

Ul (r) = uf (r) + e 2iS ^ k) u [ p (r) , (24) 

and obeys t^(0) = 0. Thus we obtain 

8i{k) = 2 Re p e (k, 0). (25) 

Furthermore, obeys a simple, non-linear differential equation obtained by substituting 
u { p into eq. ([22|), 

- iP'; - 2ikp t (kr)(3' l + 2(f3' l ) 2 + ~<7 X (r) = , (26) 

where primes denote differentiation with respect to r, and 

Pi(x) = — In (27) 

is a simple rational function of x. 

We solve eq. Q213D numerically, integrating from r = oo to r = with /^(&, oo) = 
Pi(k, oo) = 0, to get the exact phase shifts. To get the Born approximation to we 
solve the equation iteratively, writing = gfim + g 2 /3n + • • •, where /3n satisfies 

- 0l % - 2ikp e (kr)/3' ei + \ X {r) = (28) 

and satisfies 

- z/% - 2ik Pl (kr)f3' i2 + 2{(3' n f = 0.. (29) 

We can solve efficiently for and /3^ 2 simultaneously by combining these two equations 
into a coupled differential equation for the vector (Pn, fti-i). This method is much faster 

10 



than calculating the Born terms directly as iterated integrals in r and will generalize easily 
to a theory requiring higher-order counterterms. 

Having found the phase shifts, we then use Levinson's theorem to count bound states. 
We then find the energies of these bound states by using a shooting method to solve the 
corresponding eigenvalue equation. We use the effective range approximation || to calculate 
the phase shift and bound state energy near the threshold for forming an s-wave bound state. 

IV. RESULTS 

For the model at hand, we calculated the energy £[h] for a two-parameter (d and w) 
family of gaussian backgrounds 

h(r) = -dve~ r2v2/2w2 . (30) 

In Fig. |2|, we show results which are representative of our findings in general. We plot the 
energy of configurations with fixed d = 1 as a function of w, for g = 1, 2, 4, 8 (from top to 
bottom). We note that to this order, for g = 8 the vacuum is unstable to the formation of 
large = regions. 
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FIG. 2. in units of v, for d = 1 and g = 1,2, 4, 8, as a function of w. 

To explore whether the charged scalar forms a non-topological soliton in a given <fi back- 
ground, we add to £[h] the energy of a "valence" ip particle in the lowest bound state. We 
then compare this total energy to M, the energy of the ip particle in a flat background, to 
see if the soliton is favored. This is the scalar model analogue of t-quark bag formation. 

For fixed g and m, we varied h looking for bound states with energy E such that E + 
£[h] < M. However, for those values of g and m where we did find such solutions, we always 
found that by increasing w, we could make £[h] < 0, so that the vacuum is unstable, as we 
pointed out above in the case of g = 8 in Fig. ^|. Thus we find that if we stay in the g,m 
parameter region where the vacuum is stable, the minimum is at h — 0, so there are no 
nontopological solitons. 

Although we did not find a non-trivial solution at one-loop order in this simple model, 
our calculation demonstrates the practicality of our method. We can effectively characterize 
and search the space of field configurations, h(r), while holding the renormalized parameters 
of the theory fixed. The same methods can be used to study solitons in theories with richer 
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structure. 



V. DERIVATIVE EXPANSION 

Our results are exact to one-loop order. The derivative expansion, which is often applied 
to problems of this sort, should be accurate for slowly varying h(r). We found it useful to 
compare our results with the derivative expansion for two reasons: first, we can determine the 
range of validity (in d and w) of the derivative expansion; and second, where the derivative 
expansion is expected to be valid, it provides a check on the accuracy of our numerical work 
and C++ programming. Where expected, the two calculations did agree to the precision 
we specified (1 %). 

In our model, the first two terms in the derivative expansion of the one-loop effective 
Lagrangian can be calculated to be 



2„,4 



Ci = Cd + az + (3z 2 + 



9 v 

327T 2 



3 

'l + z) 2 \n(l + z)-z--z 2 



487T' s 'zr 1 + z 



where z = gxjM 2 = {h? + 2hv)/v 2 , a and (3 are cutoff-dependent constants, and £ ct is the 
same counterterm Lagrangian as we used in Sec. 2. For <f) A scalar field theory a similar result 
was first derived in [|Hj. The last term above is proportional to (dh) 2 , and is completely 
cancelled by a finite counterterm that implements the renormalization prescription of Sec. 2. 
In this prescription, counterterms also cancel the az and f3z 2 terms above. Thus the 0(p 2 ) 
derivative expansion for the effective energy, to be compared with the phase shift expression 
for £[h], is 



£Mh\ = j d 3 x -iyhf + ^(h 2 + 2hvf + 



l + z) 2 ln{l + z)-z--z 2 



■ (32) 



32tt 2 

The results of the comparison with the phase shift method can be seen in Fig. |3| for 
d = 0.25, and g = 4. A similar pattern holds in general for other values of both d and g. 
As the width becomes larger, the two results merge. This is as expected, since it is for large 
widths, and thus small gradients, that we expect the derivative expansion to yield accurate 
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results. As the width tends to zero, both results tend to zero, and the fact that the plot 
tends to 1 simply indicates that the derivative expansion result goes to zero faster than the 
phase shift result. 




FIG. 3. (£ — £de)/£ for d = 0.25, g = 4, as a function of w. 



VI. CONCLUSIONS 

We have presented a numerically tractable method for evaluating the one-loop effective 
energy of a static background field configuration in a renormalizable quantum field theory. 
Since we rely on calculating phase shifts, the method is only suitable for rotationally invariant 
(or generalized rotationally invariant) backgrounds. The model explored in the present work 
is particularly simple and does not support any solitons. However, our methods could just as 
well be applied to any renormalizable field theory. They can be used to study the one-loop 
quantum stability of field configurations in the standard electroweak model as well as various 
unified models that support monopoles, strings and the like. Topologically non-trivial field 
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configurations with maximal symmetry, like the "hedgehog solutions" in chiral models can 
be studied in this fashion. Ultimately we hope to be able to reliably determine whether 
large Yukawa couplings may yield solitons in the standard electroweak model. 
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